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ABSTRACT 

Building upon our previous MHD simulation study of magnetic channeling in radia- 
tively driven stellar winds, we examine here the additional dynamical effects of stellar 
rotation in the (still) 2-D axisymmctric case of an aligned dipolc surface field. In 
addition to the magnetic confinement parameter 77* introduced in Paper I, we char- 
acterize the stellar rotation in terms of a parameter W = V ro t/V OI b (the ratio of the 
equatorial surface rotation speed to orbital speed) , examining specifically models with 
moderately strong rotation W = 0.25 and 0.5, and comparing these to analogous 

non-rotating cases. Defining the associated Alfven radius Ra ~ T]l R* and Kepler 
corotation radius Rk ~ W^ 2 ^R^, we find rotation effects are weak for models with 
Ra < Rk, but can be substantial and even dominant for models with Ra ^ Rk- In 
particular, by extending our simulations to magnetic confinement parameters (up to 
r)* = 1000) that are well above those (77* = 10) considered in Paper I, we are able 
to study cases with Ra 3> Rk', we find that these do indeed show clear formation 
of the rigid-body disk predicted in previous analytic models, with however a rather 
complex, dynamic behavior characterized by both episodes of downward infall and 
outward breakout that limit the buildup of disk mass. Overall, the results provide an 
intriguing glimpse into the complex interplay between rotation and magnetic confine- 
ment, and form the basis for a full MHD description of the rigid-body disks expected 
in strongly magnetic Bp stars like a Ori E. 
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Stars: rotation — Stars: mass loss 



1 INTRODUCTION 



In Paper I of this series (ud-Doula & Owocki 20021, we ex- 



amined the effect of a large-scale dipole magnetic field on 
the radiatively driven wind from a hot, massive star. The 
radiative envelopes of such hot stars means they lack the 
strong convection zone that drives the dynamo generation of 
magnetic activity cycles in the Sun and other relatively cool 
stars. Nonetheless, in recent years spectro-polarimetric ob- 
servations have led to positive detections of large-scale fields 
in several such hot stars, often well fit by a dipole tilted rel- 
ative to the star's rotation axis (e.g., Donati et al.|2002 |. In 
some cases the associated period of rotational modulation is 
quite long, weeks or even years (e.g. in 8 Ori C, HD191612; 
Donati et al. 20061), implying that the direct dynamical effect 



of rotation on the magnetic channeling of the wind is likely 
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to be limited. As a first approximation, the MHD simulation 
models of Paper I thus ignored the effects of rotation. 

More generally, however, massive stars tend to have 
quite rapid rotation, as evidenced both by the substan- 



tial broadening in photospheric spectral lines (Conti & 
Ebbets 1977 Fukuda 1982[ ), which indicate projected ro- 



tation speeds of hundreds of km/s, and by the relatively 
short-period of observed modulations for some stars, e.g. 
the magnetic Bp star a Ori E, for which the inferred ro- 



tation period is about 1.2 d ( Walborn 1981 1. Both lines of 



evidence suggest that hot-star rotation rates are commonly 
a substantial fraction of the 'critical' rate at which the equa- 
torial surface would be in Keplerian orbit. Since this implies 
centrifugal forces that are comparable to the inward pull of 
gravity, it is clear that such levels of rotation could signifi- 
cantly influence the magnetic channeling of a stellar wind. 

Previous studies have focused in particular on the po- 
tential role of magnetic fields in spinning up the wind outflow 
and channeling it into an equatorial disk that might be cen- 
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trifugally supported against gravity. Cassinelli et al. ( 2002 1 



argued that such magnetic spin-up could effectively eject 
material into a "Magnetically Torqued Disk" (MTD), in 
which individual fluid elements would be in local Keplerian 
orbit, and specifically proposed this as a model for the Kep- 
lerian "decretion disks" inferred for Be stars. For the chem- 
ically peculiar Bp stars that have been directl y observed to 
have very strong magnetic fields (> 10 4 



G), 



Townsend & 



Owocki (2005 herefter TO-05) developed a somewhat dif- 



ferent "Rigidly Rotating Magnetosphere" (RRM) paradigm 
in which the field again spins up and channels wind mate- 
rial into a disk, but now is also sufficiently strong to hold 
it in rigid body rotation. This RRM model has proven par- 
ticularly successful in explaining the rotationally modulated 



Balmer emission observed from a Ori E (Townsend. Owocki 
fc Groote|2005 l. 



To test these semi-analytic paradigms, we have made 
some initial efforts to extend the numerical MHD simula- 
tions of Paper I to include rotation in the simple 2D axisym- 
metric case of a rotation-aligned dipole. The results indicate 
that a large-scale field strong enough to torque wind mate- 
rial to Keplerian orbital speed tends also to propel material 
away from the star, rather than into the kind of stationary, 



Keplerian disk envisioned in the MTD model (ud-Doula, 



Townsend fc Owocki 2006. Owoc kI7Townsend fc ud-D oula 
2006a|b|). However, for cases with strong-enough magnetic 



confinement to hold material down against such outward es- 
cape, there can indeed form a limited rigid-body disk quite 



similar to that predicted by the RRM model ( Owocki & ud- 
Doula 2003] |OwockT||2006| ). But even in such cases, there is 
irregular breakout of material from the outer disk, leading to 
a sudden magnetic reconnection heating that could explain 
the hard X-ray flares seen from a Ori E (ud-Doula et al 
2006). 



The present paper further extends these previous MHD 
simulations with a more extensive 2-D parameter study cov- 
ering models over a range in both rotation rate and degree 
of magnetic confinement. A particular focus is to develop a 
clearer physical picture of the complex competition between 
wind-fed build-up of material in a disk vs. losses by both 
outward ejection and infall back to the starQ Moreover, the 
broad parameter study here allows us to examine in detail 
how these processes are affected by various combinations of 
rotation rate and magnetic field strength. To lay the basis 
for the results presented in section 4, section 2 first reviews 
the general numerical MHD approach, and section 3 defines 
the overall parameter domain. Section 5 concludes with a 
summary and outline for future work. 



2 NUMERICAL METHOD 

2.1 Vector form of basic MHD equations 

As in Paper I, our general approach is to use the ZEUS-3D 
(Stone & Norman 19921 numerical MHD code to evolve a 



consistent dynamical solution for a line-driven stellar wind 
from a star with a dipole surface field. Our implementation 



here again adopts spherical polar coordinates with radius r, 
co-latitude 8, and azimuth <j>, but now in a "2.5-D" formula- 
tion that allows for non-zero azimuthal components of both 
the magnetic field B<p and velocity v$, while still assuming 
all quantities are constant in the azimuthal coordinate angle 
(f>. To maintain this 2.5-D axisymmetry, we assume the stel- 
lar magnetic field to be a pure-dipole with polar axis aligned 
with the rotation axis of the star. 

In vector form, the standard formulation of magneto- 
hydrodynamics includes equations for mass continuity, 

dp 



+ v ■ Vp + pV • v = , 



(1) 



and momentum balance, 
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where the notation follows common conventions, and is de- 
fined in detail in section 2 of Paper I. (Note that eqn. |2| 
here corrects some minor errors in the corresponding eqn. 
(2) of Paper I.) 



2.2 Rotation terms in the advective acceleration 

The inclusion of a finite rotation in the present work leads 
to additional nonzero terms proportional to the azimuthal 
velocity within the three vector components of the ad- 
vective acceleration, 
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(4) 



Or r 08 i i 

The terms here without derivatives represent the inertial 
forces arising from the curvature of the coordinate system, 
namely centrifugal and coriolis forces, which, e.g., in the 
absence of external torques, enforce conservation of angular 
momentum within the rotating flow described in spherical 
coordinates. 

Of course, in magnetic models the Lorentz force, rep- 
resented by the second term on the right side of eqn. Q, 
can significantly channel the flow, competing against these 
inertial terms. In rotating models, this Lorentz force can 
now also impart a significant torque to spin up the outflow; 
moreover it now includes additional terms proportional to 
the non-zero B^, which themselves represent a component 
for outward angular momentum transport. 

This competition between the magnetic Lorentz force 
and the inertia terms associated with rotation represents a 
central focus of the present study. In some ways, it paral- 
lels the central competition examined in Paper 1, namely 
between magnetic forces and the inertia associated with the 
radial outflow of the wind. 



1 To allow focus on these issues of disk buildup, we defer here 
any discussion of wind angular momentum loss and the resulting 
stellar spindown to a future, follow-up paper. 



2.3 Radial driving of wind outflow 

This radial outflow arises from the strong radial driving of 
the line-force, gimes- As in Paper I, we model this here in 
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terms of the standard Castor, Abbott & Klein ( 1975 here- 



after CAK) formalism, corrected for the finite cone angle of 
the star, using a spherical expansion approximation for the 
local flow gradients ( Pauldrach, Puis, Hummer & Kudritzki 



1985 Friend & Abbott 1986 1 and ignoring non-radial line- 



force components that can arise in a non-spherical outflow. 
Although such non-radial terms are typically only a few per- 
cent of the radial force, in non-magnetic models of rotating 
winds, they act without much competition in the lateral 
force balance, and so can have surprisingly strong effects on 



the wind channeling and rotation (Owocki, Cranmer & Gay 



Ieylp96l |Gayley fc OwocldpOOOl ). But in magnetic models 
with an already strong component of non-radial force, such 
terms are not very significant, and since their full inclusion 
substantially complicates both the numerical computation 
and the analysis of simulation results, we have elected to de- 
fer further consideration of such non-radial line-force terms 
to future studies. 

By limiting our study to moderately fast rotation, half 
or less of the critical rate, we are also able to neglect the 
effects of stellar oblateness and gravity darkening. 



2.4 Isothermal flow approximation 

Another simplification retained from Paper 1 is that the flow 
is strictly isothermal. This is generally a reasonable approx- 
imation in steady-state, spherical wind models, wherein the 
competition between photoionization heating and radiative 
cooling keeps the wind close to the stellar effective temper- 
ature (Pauldrach 1987 Drewj 1989). However, models with 
significant magnetic channeling can guide the flow toward 
strong shock compressions that heat the gas to tempera- 
tures of millions of Kelvin. The associated extensive X-ray 
emission fBabel & Montmcrlc 1997a b) has indeed been a 
focus of some our previous simulations aimed at modeling 
the observed X-ray spectrum from 9 1 Ori C (Gagne et al 



2005). Moreover, our other simulations show that reconnec- 



tion heating associated with centrifugally driven breakout 
events might provide a basis for explaining the relatively 
hard X-ray flare events seen in the magnetic B-star a Ori E 
( ud-Doula et al.|2006 l. 

While such detailed treatments of the wind energy bal- 
ance can thus be quite important for modeling the X-ray 
emission from specific stars, including this in the general 
parameter study here would require introducing an addi- 
tional free parameter, associated with the wind density, and 
representing the relative importance of radiative cooling in 
the post-shock region ( |ud-Doula|2'00 3). This would in effect 
require a i/iree-dimensional parameter study, representing 
cooling, magnetic confinement, and rotation. To maintain 
the focus here on just the one additional degree of freedom 
associated with rotation, built upon the study of isothermal 
magnetic confinement in Paper I, we again assume a simple 
isothermal wind with the wind temperature kept equal to 
the stellar effective temperature, here taken to be 50,000 K. 

A further advantage is that, at such a temperature, the 
gas pressure terms in eqn. |2| are typically unimportant 
throughout most of the supersonic outflow. Thus, although 
these pressure terms are still fully included in the numerical 
simulations, they can be largely ignored in interpretation 
of results, allowing for a focus on the dominant competing 



forces associated with gravity, radiative driving, magnetic 
field and flow inertia. 



2.5 Boundary conditions and numerical 
considerations 

Finally, numerical specifications such as the computational 
grid and boundary conditions are again similar to Paper 
I, except of course that the lower boundary now has a non- 
zero azimuthal speed v$(R*, 9) — V I0 t sin#, where Kot is the 
equatorial surface rotation speed. Also, instead of setting the 
azimuthal field to zero at the stellar surface, as in Paper I, we 
now compute a generally non-zero B$ at the lower boundary 
ghost-zone by linear extrapolation from the values in the 
two innermost zones of the actual computational grid. This 
assumes vanishing second derivatives of the azimuthal field 
components. 

In all simulations presented here the time step is based 
on relatively low Courant number of 0.3, a choice that helps 
ensure stability and reduced error in computed shock prop- 
erties (Falle 2002). Simulations of selected models with an 
even lower Courant number of 0.1 gave very similar results 
to the standard runs. 



3 TWO-PARAMETER STUDY 

3.1 Magnetic confinement parameter 77* 

Let us now consider how best to frame our parameter study 
for the combined effects of rotation and magnetic channeling 
in a line-driven wind. In the absence of significant rotation or 
magnetic fields, the line- force overcomes the stellar gravity 
to drive a nearly radial wind outflow characterized by a mass 
loss rate M and terminal wind speed V^, . When a magnetic 
field is added, the inertia of this radial outflow competes 
against the Lorentz forces. A key result of Paper I is that 
the overall net effect of a magnetic field in diverting such 
a wind outflow can be characterized by a single magnetic 
confinement parameter, 



MVoo 



(6) 



where B eq is the surface field strength at the magnetic equa- 
tor. This sets the scale of the ratio of magnetic energy to 
wind kinetic energy, 

ir/R,) 2 ~ 2q 



r){r) 



B 2 /8tv 
pv 2 /2 



= V* 



V A 



= Afi 



(7) 

The last two equalities emphasize this energy ratio can also 
be cast as the square of the ratio of the Alfven speed, Va = 
B/^/4np, to flow speed, v, i.e. as the inverse square of the 
Alfvenic Mach number, Ma = v/Va- 

The square bracket factor in the middle equality shows 
the overall radial variation; q is the power-law exponent for 
radial decline of the assumed stellar field, e.g. q = 3 for a 
pure dipole, and j3 is the velocity-law index, with typically 
/3 « 1. For a star with a non-zero field, we have n* > 0, and 
so given the vanishing of the flow speed at the atmospheric 
wind base, this energy ratio always starts as a large number 
near the stellar surface, n(r — * R,) — *■ 00. But from there 
outward it declines quite steeply, asymptotically as r~ 4 for a 
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dipole, crossing unity at the Alfven radius denned implicitly 
by n{R A ) = 1. 

For a canonical /3 = 1 wind velocity law, explicit solu- 
tion for Ra along the magnetic equator requires finding the 
appropriate root of 

2q-2 / u \ 2q-3 

= »?. , (8) 



Ra 
R, 



Ra 



which for integer 2q is just a simple polynomial, specifically 
a quadratic, cubic, or quartic for q = 2, 2.5, or 3. Even 
for non-integer values of 2q, the relevant solutions can be 
approximated (via numerical fitting) to within a few percent 
by the simple general expression, 

Ra 
R, 



l + (r ? » + l/4) 1/(29 - 2) -(l/4) 1/(29 - 2) 



(9) 



For weak confinement, rj* <C 1, we find Ra — * R*, while for 
strong confinement, 7/* ^> 1, we obtain -Ra — > r/t^ 2q ~ 2 ^ R„. 
In particular, for the standard dipole case with q — 3, we 
expect the strong-confinement scaling Ra/R* ~ 0.3 + rji^ . 

Clearly .Ra represents the radius at which the wind 
speed v exceeds the local Alfven speed Va- But Paper I 
showed that it also characterizes the maximum radius where 
the magnetic field still dominates over the wind, and is just 
somewhat above (i.e., by 20-30%) the maximum extent of 
closed loops in the magnetosphere. Moreover, as we shall 
see below, in rotating winds these closed loop regions tend 
to co-rotate nearly rigidly with the underlying star, and in 
this sense -Ra is just above the maximum radius for wind 
co-rotation near the equator. For convenience in discussing 
results, let us thus denote the maximum radius of such closed 
(and generally co-rotating) loops as 



R c ~ R* + 0.7(R A - R*) . 



(10) 



3.2 Rotation parameter W 



Let us next seek a similarly convenient parameterization for 
the stellar rotation. This can again be characterized in terms 
of a speed, namely the equatorial surface rotation speed Kot- 
But instead of relating that to the flow speed or Alfven speed 
in the stellar wind, the stellar origin of rotation suggests it 
may be better to compare it to a speed representative of the 
gravity at the stellar surface. Specifically, let us thus define 
our dimensionless rotation parameter as 



W ■ 



Vrot 
Vorb 



(11) 



where Vovb = \J GM/R, is the orbital speed near the equato- 
rial surfac^] This characterizes the azimuthal speed needed 
for the outward centrifugal forces to balance the stellar sur- 
face gravity. It is only a factor l/v2 less than the speed 



2 This is closely related to the commonly used rotation parame- 
ter uj = Q/Q cr i t , defined by the star's angular rotating frequency 
Q relative to the value this would have as the star approaches 
"critical" rotation, Q c . Our choice here more directly relates to 
the additional local speed needed to propel material into Keple- 
rian orbit, and avoids some subtle assumptions (e.g. rigid-body 
rotation using a Roche potential for gravity) about how the global 
stellar envelope structure adjusts to approaching the critical ro- 
tation limit. 



Vcsc needed to fully escape the star's surface gravitational 
potential. 

For a non-magnetic rotating star, conservation of angu- 
lar momentum in a wind outflow causes the azimuthal speed 
near the equator to decline outward as ~ 1/r, meaning 
that rotation effects tend to be of diminishing importance 
in the outer wind. 

By contrast, in a rotating star with a sufficiently strong 
magnetic field, magnetic torques on the wind can spin it 
up; for some region near the star, i.e., up to about the maxi- 
mum loop closure radius R c , they can even maintain a nearly 
rigid-body rotation, for which the azimuthal speed now in- 
creases outward in proportion to the radius, 



v<j,(r) = V TC 



R* 



<R 



(12) 



As such, even for a star with surface rotation below the or- 
bital speed, W < 1, maintaining rigid rotation will eventu- 
ally lead to a balance between the outward centrifugal force 
from rotation and the inward force of gravity, 



Rk 



GM 

R K 



(13) 



Combining this with eqns. (Ill and (12 1 gives a simple ex- 



pression for the associated "Kepler radius" 
-Rk = W" 2/3 R» . 



(14) 



Unsupported material at radii r < Rk will tend to fall 
back toward the star, but any material maintained in rigid- 
rotation to radii r > Rk will have a centrifugal force that 
exceeds gravity, and so will tend to be propelled further out- 
ward. Indeed, any corotating material above an "escape ra- 
dius", which is only slightly beyond the Kepler radius, 

Re = 2 1/3 R k , (15) 

will have sufficient rotational energy to escape altogether the 
local gravitational potential, unless, of course, temporarily 
held down by the magnetic field. 



3.3 2-D parameter grid of models 

The circles in figure [T] lay out the 2-D grid of models com- 
puted for the present study, plotted in the plane of rotation 
parameter W vs. log of the magnetic confinement parameter 
77* . The models along the x-axis include several specific cases 
already examined in the non-rotating study in Paper I, with 
now however some additional extensions toward the strong 
confinement limit, viz. log 77* = 1.5, 2, 2.5 and 3. In addition, 
there are now two new sets of corresponding models with ro- 
tation parameters W = 0.25 and 0.5. As noted, we do not 
consider faster rotation than W = 0.5 because this would 
introduce a significant stellar oblateness that would compli- 
cate specification of the lower boundary condition for the 
spherical coordinate system used in the ZEUS MHD code. 

The solid curve in figure [l] represents the parameter 
combination for which Ra = -Rk in the dipole case (q = 3) 
with velocity index f3 = 1. This contour thus roughly divides 
the parameter space diagonally: models below and to the left 
have only slow rotation and/or weak confinement, and so 
Ra < -Rk; models above and to the right have fast rotation 
and/or strong confinement, and so -Ra > -Rk. 

Our analysis of the associated simulations show that the 
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lower-left models give generally quite similar overall struc- 
ture to what was found for the non-rotating models in Pa- 
per I. The more interesting cases are those in the regions 
above and/or to the right, and in the transition region with 
Ra « Rk- 

The transition region represents cases for which the 
magnetic spin-up is just adequate to propel material into 
Keplerian orbit. As such, it might seem to be appropriately 
fine-tuned to produce the kind of magnetically torqued disk 
advocated by Cassinelli et al. (20021. However, as discussed 
below and in previous papers ( ud-Doula et al.|2006 Owocki 



ir w 



2006| pwocki et al.||20~06a|b| ), our simulations indicate that 



even for these optimal parameter cases, the rotating magne- 
tophere is characterized by a combination of infall and out- 
flow respectively below and above the Kepler radius, with no 
apparent tendency to form an extended, stable, Keplerian 
disk. 

On the other hand, in the limit of strong confinement 
with Ra S> Rk, the dominance of the field can confine the 
material in a rigid-body disk, as postulated in the "Rigidly 
Rotating Magnetosphere" (RRM) formalism developed by 
TO-05. The full MHD simulations here allow us to directly 
test this RRM concept, and define its limitations as wind 
material accumulates in the disk, leading eventually to a 
centrifugally driven breakout overcoming the confining mag- 



netic tension (see ud-Doula et al. 2006 and Appendix of TO- 



05). Toward this goal, the magnetic confinement parameters 
considered here extend to values (r/ t — 1000) that signif- 
icantly exceed the maximum (77* = 10) attempted in the 
non-rotating study of Paper I. 

Such models with strong magnetic confinement are, in 
fact, significantly more computationally challenging, since 
the greater rigidity of the magnetic field implies a higher 
Alfven speed (see Appendix), and thus requires a smaller 
numerical time step to maintain stability under the Courant 
criterion. Indeed, this problem is often exacerbated by the 
tendency for the strong, nearly horizontal field near the mag- 
netic equator to completely inhibit any wind base outflow 
there; this leads then to short-lived nearly evacuated regions 
where the Alfven speed can become exceedingly large, for- 
mally even approaching the speed of light! To keep the time- 
step from becoming too small, we thus choose to artificially 
add mass to these small evacuated regions at a level that is 
sufficient to limit the local Alfven speed to 



where 



Va ^ max(20,000 km/s, max(V Ap )) . 



max(VAp) = 0.65,/r^Kx 



(16) 



(17) 



is the expected maximum polar Alfven speed, as given by 
the analysis in the Appendix. We check that the amount of 
artificially added mass is still quite insignificant compared 
to the global mass loss in the wind, i.e. less than a percent 
in even the most extreme (77* = 10 3 ) cases. 

Moreover, the larger Alfven radius means such mod- 
els need generally a larger outer boundary radius, and the 
larger breakout timescale (as predicted in the Appendix of 
TO-05) means that models have to be run longer to cover 
the breakout cycles and associated accumulation of mass 
in any RRM disk. Finally, as discussed further below, the 
closed magnetic topology of episodic outbursts can compli- 
cate the proper specification of the outer boundary condi- 
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Figure 1. Plot of rotation vs. confinement parameter, W vs. 
log 77*, to define the models in our 2D parameter study, repre- 
sented here as filled and open circles, with the filled circles repre- 
senting models of principal focus for more detailed analyses. The 
solid curve respresents the contour for models with R^ = Rk 
when the magnetic index is set to the dipole value q = 3. Models 
above this have R\ > i?x an d thus strong confinement plus rapid 
rotation, whereas models below have _Ra < Rk and thus either 
weak confinement or relatively slow rotation. 



tion, and in practice reflection effects as these outbursts are 
advected through the boundary can occasionally even halt 
the computation altogether. In summary, the extension of 
MHD simulations into the very strong confinement domain 
remains an ongoing challenge. 



3.4 Stellar and wind parameters 

Much of the procedures in the current study follows Paper 
I. Specifically, we use the same standard non-magnetic and 
non-rotating wind model used in Paper 1, but now at the ini- 
tial time we suddenly introduce both a dipole magnetic field, 
and a surface rotation at the lower boundary, with both de- 
fined relative to a common polar axis. This standard model 
has stellar parameters representative of a typical OB super- 
giant, with a radius _R* = 19i?0, a luminosity L — 10 6 I/q, 
and an effective mass of M — 25M . (This reflects a factor 
two reduction below the Newtonian mass to account for the 
outward force from the electron scattering continuum.) This 
radius and mass imply an effective equatorial surface orbital 
speed of V OT b ~ 500 km/s. 

The line-driving assumes a CAK power index a — 0.6 
and a line normalization such that the non-magnetic, non- 
rotating wind has a mass loss rate of about M « 3 x 
10 -6 M©/yr and a terminal speed of about Voo ~ 2400 km/s. 
Since both the stellar and wind parameters are fixed, we vary 
the magnetic confinement parameter 77* solely through the 
variations in the assumed equatorial surface field strength, 
B eq . As noted above, we do not consider rotation parame- 
ters W > 1/2 since this would deform the stellar surface and 
require consideration of gravity darkening, neither of which 
are taken into account in our models. 

One further difference compared to the simulations in 
Paper I is that we find it necessary to run the rotating mod- 
els here for a longer time in order to identify properties of a 
relaxed, quasi-stationary asymptotic state (especially for the 
strong confinement models). To facilitate comparison among 
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different cases, we standardize a run to duration of t = to 
3 Msec, which is already 6 times longer than the 0.5 Msec 
used in Paper I. But we have also run selected models for a 
longer time, e.g. 6 Msec for our standard case with W — 1/2 
and 77* = 100. Required run times per model are typically 
about one to two weeks on a standard workstation. 



4 RESULTS 

4.1 Co-rotation and rigid disk 

As a standard example to frame the overall study here, let 
us first focus on this case with confinement n, — 100 and 
rotation W = 1/2 (K-ot = 250 km/s). Figures [5] show a se- 
ries of time snapshots of the 2D spatial configuration of the 
magnetic field (solid lines), with the color scale representing 
logarithm of density; figure [3] give a similar time sequence 
for the azimuthal flow speed, scaled relative to the value 
that would occur in rigid-rotation, i.e. \ = v<p(r)/£lr sin#, 
where = V I0 t/R* is the star's angular rotation frequency. 
The time snapshots were chosen to illustrate both rela- 
tively quiescent intervals (top panels), and phases with dy- 
namic centrifugal breakout (bottom panels). The dashed cir- 
cle represents the Kepler co-rotation radius at the equator 
(Rk « 1.6i?»). 

In the evolution immediately following the initial condi- 
tion, the magnetic field channels wind material toward the 
tops of closed loops near the equator, where the collision 
with the opposite stream leads to a dense disk-like structure 
(see top panels). But the gas is also generally torqued by 
the field, with, as can be seen in the upper panels of figure 
[3] material in the closed magnetosphere up to R c w 2.7 7i* 
kept nearly in rigid-body co-rotation with the star. Note that 
these closed, rigidly rotating loops thus extend through and 
beyond the Kepler radius. For any material trapped on loops 
below Rk, the outward centrifugal support is less than the 
inward pull of gravity; since much of this material is com- 
pressed into clumps that are too dense to be significantly 
line-driven, it thus eventually falls back to the star follow- 
ing complex patterns along the closed-field loops. 

By contrast, the dense material above the dashed line 
at Rk has a net radially outward force from the centrifugal 
acceleration vs. gravity. Still, during the initial build-up of 
this material at the tops of loops above Rk, the magnetic 
field provides tension that is strong enough to hold it down, 
forming then a segment of the rigidly rotating disk predicted 
in the analytic RRM analysis by TO-05. However, much 
as anticipated in the Appendix of their paper, eventually 
material in the outer region of this RRM accumulates to 
sufficient density to force open the magnetic field, leading to 
the kind of centrifugally driven breakout events simulated in 
( ud-Doula et al.|2006[ ). This is illustrated here in the bottom 
panels of figure [2| 



Note however from figure[3]that certain regions, marked 
in blue, actually have a net azimuthal motion that is against 
the sense of the stellar rotation. This surprising and coun- 
terintuitive result is not a numerical artifact, but rather is 
related to a reverse torque effect that occurs in regions of 
rapid wind acceleration. This was first discussed b y |Mac- 
Gregor fc Friend] |l987), who extended the classic |Weber 



the solar wind to the case of the more rapidly accelerating 
line-driven winds. The 2-D analog for the dipole field here 
has little impact on our study of rotation and magnetic con- 
finement, and so we defer further discussion to an upcoming 
paper that focuses on the role of the magnetic field in out- 
ward angular momentum transport and stellar spindown. 

The centrifugal breakouts occuring in these simulations 
necessarily imply a breakdown in the basic formulation for 
ideal MHD within the Zeus code. In equatorial regions where 
the wind or centrifugal terms stretches out field lines of op- 
posite polarity, the finite grid resolution allows effective re- 
connection, with the associated release of magnetic energy 
effectively lost instantaneously (e.g. due to radiative cooling) 
in these isothermal simulations. This is admittedly a simpli- 
fied representation of the very complex physics thought to 
occur in actual reconnection, which indeed is an intense area 
of modern plasma physics research ( Shay et al.||1999| ). But 
in the present context of driven reconnection, the overall 
global evolution seems likely to be set by the central compe- 
tition between magnetic confinement and centrifugal break- 
out, with relatively little sensitivity to the details of local 
reconnection sites. 



4.2 Global evolution of equatorial disk in radius 
and time 

A key result of the simulations here is that there is really 
no true steady state possible, since the secular buildup of 
material in the disk must eventually lead to an episodic ma- 
terial breakout once the centrifugal forces overwhelm the 
finite magnetic tension. One primary goal of the more ex- 
tensive parameter study here is to examine in detail the 
nature of this build-up and dissipation of mass in an RRM 
disk, and how this varies with the changes in the rotation 
rate and magnetic confinement. To facilitate illustration of 
these competing processes, let us define a radial mass dis- 
tribution of the disk, computed at each radius r in terms of 
the mass within some specified co-latitude range about the 
equator, 



dm e {r,t) _ 2 ^ 2 
dr 



Tr/2+AB/2 



/2-A0/2 



p(r,e,t) sin6 d8. (18) 



& Davis ( 1967 1 1-D magnetic monopole rotation model for 



To isolate the disk but not miss too much disk material 
during various oscillations about the equator, we choose a 
narrow, but not-too-limited range AG = 10°. Figure [4] shows 
schematically how this is computed. 

For the same standard case (W = 0.5, 77* = 100) shown 
in figures [2] and [3] figure [5] shows a colorscale plot of this 
disk mass vs. radius (on the ordinate) and time (on the 
abscissa). The horizontal lines mark, from top to bottom, 
the estimated Alfven radius 7?a ~ 3.47?*, the loop closure 
radius R c , Kepler radius Rk ~ 1.67?*, and inner disk radius, 
R m = (2/3) 1/3 7? K « 1.47?* [see eqn. (19) of TO-05]. Within 
the RRM model, the last represents the location where the 
effective potential along a rigid-field loop first develops a 
local minimum, which can then trap material fed from the 
wind. 

The plot shows quite succinctly, and vividly, the global 
time evolution of the equatorial disk material. Initially mass 
builds up in the region around Rk, but then there appear 
repeated episodes of infall of inner disk material back onto 
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102 ksec 402ksec 702 ksec 




Figure 2. Snapshots of density (in cgs units on a logarithmic color-scale) and field lines (solid lines) at the labeled time intervals for the 
model with r;* = 100 and W = 1/2. The top panels show the model during a relatively quiescent period, when a dense rigidly rotating 
disk is being gradually built up. Note, however that material below the Kepler co-rotation radius shown as a dashed circle) falls 

back onto the stellar surface, due to the lack of sufficient centrifugal support. The extent of the disk in this phase is determined by the 
magnetic field strength and extends up to the Alfven radius R\ ps 3.4JJ*, which is somewhat above the maximum outer radius of closed 
magnetic loops, R c . The bottom panels show the model later in the evolution, during one of the episodic centrifugal breakout events. 



the star, about every 200 ksec or so. This leads to a gradual 
outward progression to the lower edge of the disk material. 

But over a somewhat longer timescale, about every 
1 Msec or so, there appears another, somewhat different 
kind of disruption, one that starts higher up, closer to the 
closure and Alfven radii. This is characterized by outward 
ejection of the upper disk mass, but then also a "rebound" 
that propagates back down toward the Kepler radius, push- 
ing the trapped disk material inward, and inducing a further 
leakage of disk mass through downward infall. The overall 
effect is to regulate the disk mass so that, by the end of the 
simulation at t = 3 Msec, the addition of new material from 
the wind becomes roughly balanced by the losses to both 
infall and ejection. 

The overall structure is certainly quite dynamic, but 
near the Kepler radius there nonetheless appears to be a 
quasi-permanent disk segment that corresponds roughly to 
what is predicted by the RRM analytic analysis of TO-05, 
as well as by the recent time-dependent Rigid-Field Hydro- 



dynamics (RFHD) simulations of Townsend, Owocki, and 
ud-Doula (2007). 

4.3 Comparison with RRM breakout analysis 

Both the RRM and RFHD approaches are based on the ide- 
alization that the field is arbitrarily strong, and so remains 
perfectly rigid regardless of the amount of material in either 
the wind outflow or the disk buildup. But the associated 
discussion for both approaches recognized that the secular 
accumulation of material in the rigid-body disk would even- 
tually cause the outward centrifugal forces to overwhelm the 
available inward tension associated with any large but fi- 
nite magnetic field. In fact, the Appendix of TO-05 presents 
a simplified, but quantitative analysis of the resulting ex- 
pected "breakout" of accumulated disk material. This an- 
ticipates, at least in general terms, several aspects of the 
processes seen in the present MHD simulations. In particu- 
lar, it makes quite specific predictions for both the breakout 
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Figure 3. For the same model and same time snapshots as in figure [2] the azimuthal speed scaled by the local co-rotation speed, 
X = v^/Qrsind) . During the quiescent period the magnetosphere extends nearly to the closure radius R c 2.7R*, and co-rotates almost 
rigidly with the star (top panels). However, during episodic centrifugal breakouts (bottom panels) this rigidly rotating magnetosphere 
shrinks nearly to Rk ~ 1.6B*, represented by the dashed circle. In addition to the mass lost to the breakout, this also helps some of the 
disk mass to leak inward as infall back onto the star. 



Star 




Figure 4. Schematic diagram illustrating the computation of 
the radial mass distribution of the equatorial disk, dm e /dr (see 
eqn. \18\ . We choose a cone angle of A9 = 10° centered on the 
magnetic equator to encompass most of the material in the disk. 
We then compute the total mass Am contained in a narrow strip 
Ar within the cone. Note, however, that some of the material may 
not lie within the cone during infall episodes. 



timescales as a function of radius, and for the asymptotic 
mass accumulation near the Kepler radius. In this section, 



let us thus attempt a specific, semi-quantitative compari- 
son between those predictions and the results of the MHD 
simulations. 

4-3.1 Breakout timescale 

Without benefit of the global dynamical picture available 
from the MHD simulations here, the TO-05 breakout anal- 
ysis focused instead on the conditions for breakout at each 
local radius, conveniently scaled in terms of the Kepler ra- 
dius as £ = r/i?K- From TO-05 eqn. (A6) we find that, in 
terms of the free-fall time tg = R*/V ESC , the breakout time 
for some scaled outer disk radius £ is given by 

tb ~ rj,tft 3 — . (19) 

So t 

Here we have approximated y/iv/fj,* w 2 (cf. eqn. A8 of TO- 
05), and used the ratio of wind terminal speed to escape 
speed, Voo/Vosc ~ 3, to convert the disk confinement param- 
eter (also denoted 77,) in TO-05 into the wind confinement 
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Figure 5. For MHD simulations of the standard model case with 
r\ t = 100 and W = 1/2, the logarithm of the radial distribution of 
mass, dm e /dr, within a cone of width A9 = 10° centered on the 
equator, plotted vs. time and radius (in units of stellar radius). 
The color bar shows log(dm e /dr) in units of Mq/R*. The hori- 
zontal lines indicate the Alfven radius R\ (solid), the maximum 
loop closure radius R c (dotted), the Kepler radius i?x (dashed), 
and the inner RRM disk radius R m = (2/3) 1 / 3 R K (dot-dashed). 
Note the accumulation of mass in a disk near the Kepler radius, 
and how this is limited both by relatively frequent (ca. every 
200 ksec) intervals of downward infall, and by less-frequent (ca. 
every 1 Msec) major disruptions with substantial upward break- 
out. 



parameter defined her^] Note also that the Kepler-scaled 
stellar radius can be written as = ii*/i?K = W 2 ^ 3 '. 

Applying our stellar free-fall time iff w 19 ksec, then 
for our standard (?7* = 100, W = 1/2) model, the predicted 
breakout time is 



7Msec 



(20) 



In the MHD simulations for this standard case, the break- 
outs seem to originate around R c ~ 2.7i?*, for example as 
indicated in figure [5] by the "bi- furcations" between upward 
and downward tracks that start at r < 3R, for times around 
800, 1600, and 2700 ksec. If we thus approximate the outer 
disk radius by this maximum loop closure radius, we find 
£o = Rc/Rk ~ 1.7 and so tb ~ 1 Msec, about the timescale 
between major breakout eruptions seen in these same MHD 
simulations. 

The TO-05 breakout analysis envisioned a hierarchy of 
breakout timescales, with more frequent eruptions occurring 
at larger radii; but its concluding paragraph also anticipated 
(partly based on early versions of the MHD simulations de- 
scribed here) that breakouts originating within r < 2Rk 
could also lead to substantial disruption of the entire mag- 
netosphere. The simulations here do indeed show such major 



3 Although a footnote in the Appendix of TO-05 seems to imply 
that B„ and B eq are distinct, they are in fact both equal to the 
field strength at the equatorial surface. 



Figure 6. Cumulative mass vs. time in MHD simulations for 
the standard model case. The top curve shows the total mass in 
the entire grid, i.e. including both disk and wind material. The 
middle curve shows the disk mass within 5° of the equator, and 
within radii bounded by the inner RRM disk radius Ri n and an 
outer radius R = Rc- The bottom curve shows equatorial mass 
above this outer radius. The horizontal dashed line indicates the 
cumulative disk mass predicted by the RRM breakout analysis, 
using this same outer radius. Note that it matches quite well the 
varying asymptotic value of the middle curve, showing the disk 
mass in the MHD simulation. 



disruptions, but even after these there remains substantial 
mass near the Kepler radius. 



4-3.2 Accumulated disk mass 

For a disk with a scaled outer breakout radius £ , eqn. (A10) 
of TO-05 predicts a specific scaling for the total asymptotic 
disk mass, which in terms of the parameters here can be 
written, 



m d (£ ) = 



■S^Mts T)* W 4/3 
£o 2 (£o 2 +£o + l) 



= 2.1 x 10 J M P 



£o 2 (£o 2 +Co + l) 



(21) 



where the latter gives the numerical scaling for the stellar 
and wind parameters used here. If we then apply the con- 
finement and parameters of our standard model, and use the 
characteristic breakout radius £ ~ 1.7 adopted above, we 
obtain a predicted total disk mass of mj w 1.6 x 1O _8 M0. 

Figure [6] compares this predicted mass (red horizontal 
dashed curve) with the time variation of three types of cumu- 
lative mass in the standard model MHD results. Specifically, 
the top curve (black) shows the total integrated mass in the 
entire grid, including regions of both wind and disk; the mid- 
dle curve (red) gives the cumulative mass in the equatorial 
disk region between the inner radius Ri n « 0.87/Zk ~ 1AR* 
and outer radius R a ss R c ~ 2.7R*; and finally, the bottom 
curve (blue) shows the equatorial mass above this outer disk 
radius, r > R m R c . 

Note in particular that the predicted asymptotic disk 
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mass agrees remarkably well with this MHD simulation disk 
mass as shown by the middle, red curve. Of course, this is 
partly fortuitous, since just a 10% change in the choice of 
outer disk radius £ would imply a ca. 40% change in disk 
mass. But the overall, order-magnitude agreement seems 
likely to be quite robust, and so provides a nice consistency 
check for both the RRM analysis and numerical simulations. 

Comparison of the middle and upper curves in figure [6] 
further shows that the total mass in the thin, radially limited 
disk represents about half of the total mass in the entire 
model. In part, this reflects the fact that, for such a strong 
confinement model, a large fraction of the outgoing wind 
gets channeled into the disk. But another factor is that the 
wind material is flowing outward at a very high speed, and 
so has a much shorter "residency time" than the trapped, 
relatively static material in the disk. 



4-3.3 Limitations of a localized breakout description 

Despite this general success of the local breakout analysis 
in matching both the overall breakout timescale and the 
accumulated disk mass of this MHD simulation, the specifics 
of the dynamical evolution seen in the simulation make clear 
that the breakout process is really a global phenomenon. 
As the accumulation of material in the outer disk regions 
stresses and eventually overcomes the inward restraint of the 
magnetic field, the associated outward stretching alters the 
global field, including in the inner regions near the Kepler 
radius. Moreover, once a breakout occurs, the release of this 
stretching causes the inner, closed field lines to snap back 
inward, much like a stretched rubber band after release. The 
overshoot can push disk material below the Kepler radius 
and trigger infall back onto the star. 

Overall, the wind-fed accumulation of disk mass is thus 
balanced not just by ejection outward, but also by infall in- 
ward. In contrast to the idealized picture of the breakout 
analysis, which formally predicts the timescale for breakout 
(and thus emptying) of material right at the Kepler radius 
to become arbitrarily long, the dynamic oscillation and as- 
sociated inward spillage of material limits the asymptotic 
mass accumulated in this region. This new perspective on 
the dynamical nature of the disk mass budget has potentially 
important implications for modelling and interpreting obser- 
vational diagnostics ( see|Townsend et al. 2005; Townsend, 



Owocki fc ud-Doula|2007l 



4.4 Comparison with non-rotating model 

To demonstrate further the role of rotation in how magnetic 
fields influence a wind outflow, let us now compare the re- 
sults of this 77, = 100, W = 1/2 case with the corresponding 
non-rotating model. Figure [7] illustrates the dynamic evo- 
lution of equatorial mass for this strong-confinement case 
without rotation. Comparison with figure [5] shows that there 
are still both breakout and infall episodes, but now with 
both originating from nearly the same location, at roughly 
the loop closure radius R c . This infall from throughout most 
of the closed field region reflects the lack of any centrifu- 
gal support against gravity, and as such, there is no longer 
any accumulation of material into a circumstellar disk. The 
breakouts remain, but instead of being driven by centrifugal 



3 



1„ 



1500 

Time 



Figure 7. As in figure[5] the logarithm of the radial distribution 
of equatorial mass, dm e /dr, again plotted vs. time and radius 
for a strong confinement (77, = 100) model, but now with no 
rotation (W = 0). The horizontal lines indicate the Alfvcn radius 
i?A (solid) and loop closure radius R c (dotted). Note there are 
again complex patterns of breakout and infall, but now with more 
irregular timescales, and with the infall extending up to regions of 
breakout, near the loop closure radius. As such, there is no longer 
an extended region of mass accumulation into a circumstellar disk. 



forces, these are now the result of entrainment of the field 
with the outflowing wind. The timescales for both break- 
out and infall are comparable to the rotating case, but seem 
somewhat more irregular. Overall, without the buildup in 
the disk, there is significantly less mass in the magneto- 
sphere than in the rotating case. 



4.5 Results for variation in parameters 

To provide a wider context to these detailed results for spe- 
cific cases with strong confinement, let us now analyze re- 
sults spanning a broader range of the 2-D parameter space 
for rotation and confinement strength. Before focusing fur- 
ther on equatorial disk material, let us first briefly consider 
the effects of various magnetic field strengths and rotation 
rates on the global mass loss from the stellar wind. 



4-5.1 Global mass loss 

As discussed in Paper I, one general effect of a strong mag- 
netic field on a wind is the confinement and inhibition of 
the outflow within a belt around the magnetic equator. For 
a dipole field line that reaches up to the maximum closure 
radius R c , the co-latitude 8 C at the surface footpoint satisfies 



sin6i c = s/RjRc 



(22) 



The fraction of the stellar surface area that is covered by 
closed field lines is given by cos 8 C , leaving thus only the 
remaining fraction 1 — cos 9 C as the source of wind mass 
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the upper two curves in the lower panel use a modified form 



of eqn. (23 1, 




Figure 8. Top panel: Average mass loss rate (in units of the 
non-rotating, non-magnetic case) for MHD models with various 
rotational rates W = 0, 1/4, 1/2 plotted against ^Jr\i oc B*. Note 
that for W = 1/2 the mass loss rate is nearly constant for models 
with high confinement 77, > 25. Bottom panel: Corresponding 
mass loss rates from the analytic scaling formula in eqn. pit. 



loss. For the non-rotating case, we can thus use this open- 
field fraction to estimate the overall magnetic reduction in 
the mass loss rate, 



Mb 
M b =o 



1 - v/1 - R*/Rc 



(23) 



where R c is evaluated from eqn. ( 10 1, with the Alfven radius 
Ra from eqn. |9| using q = 3. Note that this ignores higher 
order effects, such as the reduction of the mass flux in open 
field regions due to the tilt of surface field relative to the 
radial direction for wind driving (see |Qwocki fc ud-Doula| 



2004) 



The two panels in figure [8] compare the mass loss rate 
vs. ^JrjZ (oc B) for both simple analytic scalings (bottom) 
and results of numerical MHD simulations (top), with the 
lower, middle, and upper curves in each panel correspond- 
ing to the W = 0, 1/4, and 1/2 rotation models. For the 
non-rotating case, the numerical and analytic results shown 
in the lower curves are in good overall agreement. But for 
the rotating case, the upper curves in the top panel show 
that the tendency of the strong field to reduce the overall 
mass loss rate is somewhat compensated by faster rotation, 
and in the W = 1/2 case, even flattens to nearly constant 
toward the limit of strong confinement. This reflects the ad- 
ditional effect of centrifugal forces in driving the breakout 
of material initially trapped in closed loops near and be- 
low the confinement radius J? c . In effect, the rotation allows 
eventual breakout from loops that are some factor times the 
Kepler radius, say 2Rk- 

To take this into account in an analytic scaling formula, 



M b =c 



1 - y/l - R,/R c + 1 - s/l-0.5*R*/R K , (24) 



which effectively sums separate contributions from polar 
opening and rotational breakout, with the closure and Ke- 
pler radii computed from eqns. |9J, (10 1, and (14 1. With 



this generalized scaling, the overall variations of the analytic 
curves in the lower panel roughly match the corresponding 
MHD results in the upper panel. 

4-5.2 Equatorial mass and disk 

Let us next examine how the equatorial disk region is af- 
fected by variations in the magnetic confinement and ro- 
tation parameters. Figure [9] compares the radius and time 
evolution of the equatorial mass, dm e /dr, for an mosaic of 
models with various 7/* and W. The comparison provides 
a global overview of how the equatorial mass evolution is 
affected by changes in confinement and rotation. 

For weak rotation and confinement cases in the lower 
left panels, material generally escapes outward without 
much infall, with only a modest rotational enhancement in 
mass loss. But most other models again show a complex 
competition between infall and breakout, with the latter al- 
ways being less frequent and stronger. 

In particular, this complex combination of infall and 
breakout also dominates the -Ra ~ Rk transition models 
identified in figure[l] i.e. the ones here with log??* = 1/2 and 
W = 1/2 or log 77, = 3/2 and W 7 = 1/4. Such models might 
seem optimally fine-tuned to propel material into Keplerian 
orbit, and yet they show no apparent tendency for material 
to accumulate into the extended, Keplerian disk envisioned 



in the MTD scenario suggested by Cassinelli et al. (2002) 



The lack of a sharp outer cutoff in the large-scale dipole field 
makes it incompatible with the shear of a Keplerian disk, 
and without the closed loops that hold down a rigid disk in 
the strong-confinement limit, material is propelled outward 
to escape, rather than into a stable Keplerian orbit. 

As expected, accumulation into such a rigid-body disk is 
strongest for the fastest rotation, and strongest confinement, 
as shown by models at the upper right. For strong confine- 
ment but slow or no rotation, the material infall comes from 
a greater height, set by the closure radius, which increases 
roughly as R c ~ rjV A . 

This larger infall height seems also to lead to a 
somewhat longer infall timescale. Likewise, the breakout 
timescale also seems to increase for models with stronger 
confinement paramater 77*, but not quite in the linear pro- 



portion that might be suggested by eqn. ( 19 1. The reason is 



that the scaled outer confinement radius £0 appears in the 
denominator, and since this also increases with confinement 
through it dependence on the closure the radius R c , the net 
effect is to weaken the 77* sensitivity of %, especially in the 
strong confinement limit 77, ^> 1. 

Note also that the W = 1/2 rotation model with the 
strongest confinement, 7/* = 1000, is relatively stable, with- 
out the repeated equatorial infall events seen in other mod- 
els. Instead of the extensive north-south disk oscillations 
seen in other models, in this case the variations of the equa- 
torial disk are mostly symmetric about the equator, and thus 
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Figure 9. Logarithm of radial distribution of equatorial disk mass, dm e /dr vs. radius and time, for a mosaic of models with magnetic 
confinements log(r?„) = 1/2, 1, 3/2, 2, 5/2 and 3 (columns from left) and rotations W = 0, 1/4, and 1/2 (rows from bottom) with the 
same contour levels and the ranges for time and radius as in figure [5] Again, the horizontal lines indicate the Alfven radius Ra (solid) 
and the Kepler radius _Rk (dashed). 



do not induce as much "spillage" back onto the star. The 
recent analysis of "Rigid-Field Hydrodynamics" (RFHD) 



models by Townsend et al. (2007) show that both types 



of oscillation modes are allowed, with the one dominating 
in simulations depending on subtle details of the excitation 
processes. 

But overall, it seems that the basic principals gleaned 
from the detailed study of the standard, strong-confinement 
case can be quite logically generalized to understand the 
trends in properties seen from this mosaic of models span- 
ning a broad range of rotation and confinement parameters. 



5 SUMMARY & FUTURE WORK 

This paper examines the effects of field-aligned rotation on 
the magnetic channeling and confinement of a radiatively 
driven stellar wind. It builds upon the non-rotating models 
of Paper I, extending them to cases of much stronger mag- 
netic confinement (up to 77* = 1000), and comparing a full 
spectrum of models ranging from weak to strong confine- 
ment at equatorial rotation rates from zero to a substan- 
tial fraction (W =1/4 and 1/2) of the critical (or orbital) 
limit. As an initial study, it ignores the effects of oblate- 
ness, gravity darkening and limb darkening, and is based on 
an idealization of isothermal flow driven by a purely radial 
line-force. 

The key results can be summarized as follows: 

• The 2-D parameter space represented by rotation W 
and magnetic confinement rj t can be conveniently divided 



by comparing the Kepler radius Rk ~ W 2 ^R, with the 
Alfven radius R\ w ^ 1//4 _R». 

• Models with Ra < Rk have weak rotation and/or mag- 
netic confinement, with the effects of rotation limited to 
some modest enhancement in equatorial density and overall 
mass loss rate, relative to non-rotating cases. In general any 
magnetically confined material falls back to the star. 

• Transition models with _Ra ~ -Rk have a complex com- 
bination of inner region infall and outer region breakout, but 
show no signs of accumulation of material into the kind of 
extended Keplerian disk posited in the magnetically torqued 
disk paradigm of Cassinelli et al. (2002). 

• In models with _Ra > the strong magnetic confine- 
ment combines with a sufficiently rapid rotation that can 
support material against infall near and above the Kepler 
radius, with then the magnetic field both holding material 
in nearly rigid-body rotation, and keeping it confined against 
the tendency for the centrifugal force to propel material out- 
ward against gravity. 

• Such strongly confined rotation models show a clear 
accumulation into a rigid-body disk, much as predicted in the 
analytic Rigidly Rotating Magnetosphere (RRM) formalism 
of TO-05. 

• However, the present MHD simulations show that such 
rigid-body disks can be highly dynamic and variable, with 
mass accumulation regulated by a complex combination of 
inner disk infall and outer disk breakout. 

• Nonetheless, application of the simple breakout analy- 
sis introduced by TO-05 - slightly modified to identify the 
outer disk limit with the maximum radius of loop closure - 
provides a quite good, semi-quantitive agreement with both 
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the breakout timescale and limiting disk mass of the MHD 
simulations. 

• The breakout events here are similar to those discussed 
by ud-Doula, Townsend, and Owocki (2006), except the as- 
sumption here of isothermal flow precludes us from following 
that paper's specific application to modelling X-ray flares. 

There thus remains much work for future extensions of 
the present study, including relaxation of the assumptions 
of isothermal flow, purely radial driving, 2-D axisymmetry 
with aligned dipole field, and moderate rotation rates for 
which stellar oblateness can be neglected. 

To more fully test the RRM paradigm, there is also 
a need to extend MHD simulations even further into the 
large confinement limit, to approach as closely as possible 
the estimated r], w 10 7 appropriate for Bp stars like a Ori E. 

But even within the context the present set of models, 
we have ignored here another key effect of magnetic rota- 
tion, namely the outward angular momentum loss in the 
stellar wind, and the associated spindown of the underly- 
ing star. This omission was made to allow a more directed 
focus on the already quite interesting, complex, and subtle 
effects of rotation on magnetic confinement and disk forma- 
tion. But we have already carried out an extensive analysis 
of the implications of our study for angular momentum loss 
and stellar spindown, and so intend this to be the subject 
for the next paper in this series. 
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APPENDIX A: 
SPEED 



MAXIMUM POLAR ALFVEN 



We derive here eqn. (17 1 cited in section 3.3, giving a sim- 



ple scaling relation between the expected maximum polar 
Alfven speed and the wind magnetic confinement parame- 
ter. We begin by repeating eqn. |7|) relating the ratio of ma- 
gentic to wind energy density to the ratio of Alfven speed 
to flow speed, 



B 2 



pv 2 /2 



(Al) 



Here the density p can be eliminated by noting that, for a 
steady-state, magnetically channeled flow, the lack of diver- 
gence in both the field and mass flux requires, 



4npv 
B 



— constant 



M 
BpRl 



(A2) 



where we have normalized in terms of a polar surface field B p 
and global mass flux M. For a polar field that declines with 
radius as B — B p (R t /r) q , this gives for the polar Alfven 
speed, 



VX P = 



B P R 2 
M 



Bv = 477* 



R t 



(A3) 
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where the latter equality casts this in term of the magnetic 
confinement parameter, with the factor 4 correcting for the 
fact that our standard 77, is defined in terms of the equatorial 
surface field B eq , which is half the polar value B v . If we 
further assume a standard /3 — 1 velocity law, then 

/ R \ 1/2 /R \ q/2 
Va p = 2^V oc (i-^) f^J . (A4) 

Setting the radial derivative of this to zero shows that 
the maximum occurs at radius r = (1 + l/q)R,, i.e. at 
r = (4/3)i?* for the standard dipole case with q — 3. Plug- 
ging this radius into eqn. < |A4[ ), we find that the expected 
maximum Alfven speed for dipole expansion of the polar 
wind is 

max(VA P ) = ^ spq* Voo » 0.65 ^ V x ■ (A5) 

o 

This thus provides a simple rule for the maximum Alfven 
speed to allow in running MHD simulations with increasing 
magnetic confinement parameter, as discussed in section 3.3. 



